First order transition in a three dimensional disordered system 
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We present the first detailed numerical study in three dimensions of a first-order phase tran- 
sition that remains first-order in the presence of quenched disorder (specifically, the ferromag- 
netic/paramagnetic transition of the site-diluted four states Potts model). A tricritical point, which 
lies surprisingly near to the pure-system limit and is studied by means of Finite-Size Scaling, sep- 
arates the first-order and second-order parts of the critical line. This investigation has been made 
possible by a new definition of the disorder average that avoids the diverging-variance probability 
distributions that plague the standard approach. Entropy, rather than free energy, is the basic 
object in this approach that exploits a recently introduced microcanonical Monte Carlo method. 

PACS numbers: 75.40.Mg, 75.40.Cx, 75.50.Lk, 05.50.-|-q 
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The combination of phase coexistence and chemical 
disorder plays a major role, for instance, in colossal mag- 
netoresistance oxides [1]. In these situations one faces a 
fairly general question: which are the effects of quenched 
disorder 0/ on systems that undergo a first- order phase 
transition in the ideal limit of a pure sample ? For D — 3 
systems, D being the space dimension, we only know 
that disorder somehow smoothes the transition. More is 
known in D = 2, where the effects of disorder are so strong 
that the slightest concentration of impurities switches the 
transition from first-order to second-order [3, 0, 0] ■ 

An useful physical picture in Z? = 3 is provided by 
the Cardy-Jacobsen conjecture [4]. Consider a ferromag- 
netic system undergoing a first order phase transition for 
a pure sample. Let T be the temperature while p is the 
concentration of magnetic sites. A transition line, Tc{p) 
separates the ferromagnetic and the paramagnetic phases 
in the {T,p) plane. In D 3 a critical concentration is 
expected to exist, 1 > pt > 0, such that the phase tran- 
sition is of the first-order for p > pt and of the second 
order for p < pt (at pt one has a tricritical point) . When 
p approaches pt from above, the latent-heat and the sur- 
face tension vanish while the correlation-length ^(Tc(p)) 
diverges. The Universality Class is expected to be re- 
lated with that of the Random Field Ising Model (RFIM). 
However, the Cardy-Jacobsen conjecture relies on a map- 
ping between two still unsolved models (in D = 3), the 
(large Q) disordered Potts model Q and the RFIM. 

Numerical simulation is an important tool for theo- 
retical investigations in Z? = 3. In this way, large por- 
tions of the transition line Tc{p) were found to be sec- 
ond order 0, [EUl- However, the study of the tricritical 
point as well as that of the first-order part of the transi- 
tion line seemed hopeless. The problem comes from the 
long-tailed probability distribution functions (PDF) en- 
countered at Tc{p), when comparing the specific- heat or 
the magnetic susceptibility of different samples [3 . Long 
tailed PDFs follows from the standard definition of the 



quenched free-energy at temperature T as the average of 
the samples' free-energy at the same T Q, which is dom- 
inated by rare events [20|. Furthermore, the simulation 
of a sample of linear size L with previous methods is in- 
trinsically hard even for a pure system (see [13] )• In fact, 
previous work 0,11] was limited to L < 25. 

Here, we study for the first time the tricritical point 
separating the first and the second order pieces of the 
transition line. Furthermore, we characterize a first or- 
der transition that remains so in the presence of quenched 
disorder. This has been made possible by two alterna- 
tive methods of performing the sample average that avoid 
long-tailed PDFs, reproduce the correct thermodynamic 
limit, and provide complementary information. Essential 
for this study has been the capability of studying directly 
the entropy, using a recently proposed microcanonical 
Monte Carlo method combined with a cluster algo- 



rithm We studied systems of size up to L = 128, 



which allowed a neat Finite-Size Scaling investigation of 
the elusive tricritical point. 

Specifically, we consider the site diluted Q = 4 Potts 
model with periodic boundary conditions. The spins 
CTi — 1, . . . ,Q occupy the nodes of a cubic lattice with 
probability p. We consider nearest neighbor interaction: 
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The ei are quenched occupation variables, Uj = or 1 
with probability 1—p andp respectively) [2111. The pure 
system, p=l, undergoes a first order phase transition [^, 
13[ which is generally regarded as very strong. 

We introduce a real- valued conjugated momentum per 
occupied site, tt^ The total Hamiltonian is Ti. = 

WP™ + y]j eiirf /2 (the internal energy density will be e = 
H/N [ij). In the canonical ensemble, (e)T = l/(2r) + 
(WP™/iV)T. We consider instead the microcanonical en- 
semble for the extended model {ai, tt;} at fixed e, and in- 
tegrate out the {TTi} to obtain a Fluctuation-Dissipation 



2 



formalism. The basic quantity is a function of e and the 
spins, /3=(iV — 2)/(A^e — Ti'^P'") . Its microcanonical mean 
value /3|e}(e) — {13) e is the e-derivative of the entropy per 
spin, s(e), for that particular sample {e}. 

Connection with the canonical formalism is made by 
solving the equation /3{g-j(e) — 1/T = 0, that yields the 
internal energy as a function of temperature. Ther- 
modynamic stability requires j3^^^{e) to be a decreasing 
function of e. Yet, at phase coexistence and for finite 
N , it is not (see Fig. and Ref. the equation 

/3{£}(e) — 1/T = has several roots. For T = Tc, we 
name respectively ea and Bq the rightmost and leftmost 
solutions, that correspond to the energy densities of the 
coexisting disordered and ordered phases. The critical 
temperature is fixed by Maxwell construction: the e- 
integral of /3{g}(e) — l/T^ from Cd to vanishes pst . 
The surface-tension, Z', is L^^^ /2 times the integral of 
the positive part of /3{g}(e) — l/Tc for Eq < e < Cd- 

For a disordered system, one analyzes the set of func- 
tions /3{e}.(e) corresponding to a large enough number of 
samples. There are two natural possibilities. On one 
hand, one can use the Maxwell construction for each 
sample, extracting Tc, ed, Co and S and considering af- 
terwards their sample average or even their PDF, Fig. [51 
The second alternative is to compute the sample-average 
/3(e) = /3{g}(e), and then perform on it the Maxwell con- 
struction (i.e. take the sample average of s(e), rather 
than the average of the free-energy at fixed T) . 

We have empirically found that the two sample- 
averaging are equivalent in the first-order piece of the 
critical line. This is hardly surprising, because the inter- 
nal energy as a function of T is a self-averaging quan- 
tity for all temperatures but the critical one. Therefore, 
also ed, Go and Tc are self- averaging properties in the 
first-order piece of the critical line. The first method of- 
fers more information but it is computationally more de- 
manding (it requires high accuracy for each sample) . The 
method featuring /3(e) can be used as well in the second- 
order part of the critical line, nevertheless its merit in 
that region are yet to be researched. 

We have investigated the phase transition for several 
p values in the range 0.75 < p < 1. As a rule, we found 
that at fixed p the latent heat is a monotonically decreas- 
ing function of L, Fig. [3l For each p value, we simulated 
L = 16, 32, 64 and 128 (for a given p, we did not con- 
sider larger lattices once the latent heat vanished). For 
all pairs (T,p) we simulated 128 samples. Besides, some 
intermediate L values were added for the Finite Size Scal- 
ing study below (see Fig. [4]), and we have raised to 512 
the number of samples for (L = 16, 32, p = 0.86, 0.875). 

We used a Swendsen-Wang (SW) version of the micro- 
canonical cluster method [13|. For disordered systems, 
SW updates properly loosely connected regions jl6| and 
does not require painful parameter tunings. For each 
sample, we simulated at least 20 e values in the range 
— 1.2 < e < —0.5. The values of e were decreased se- 
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FIG. 1: (color online). Sample-averaged e-derivative of the 
entropy, /3(c), for several lattice sizes, L, and spins concentra- 
tions, p. Metastability requires a non-decreasing /3(e). The 
horizontal line marks the critical (inverse) temperature 1/Tc, 
obtained through Maxwell's construction. At fixed L the sur- 
face tension increases for growing p. Note that, for fixed dilu- 
tion, a seemingly first order transition [L — 64, bottom-right), 
may actually be of the second order if studied on larger lat- 
tices (L = 128, bottom- left). 
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FIG. 2: (color online). Histograms for the sample-dependent 
latent-heat A^^je = e^- eo (left) and surface-tension (right). 
In the top panels we show results in the largest lattice, where 
two very close spin concentrations behaves very differently. 
The three types of drawn horizontal lines (indicating central 
value and statistical error) correspond, from top to bottom, 
to the median, the mean and the value obtained from /3(e). 
In the lower panels we show the histograms for p = 0.98 and 
several L (mind the difference in the horizontal scales with 
the upper part). The latent- heat is self-averaging while the 
surface tension is not. 



quentially, to make use of the thermalization effort at 
the previous energy density. The microcanonical cluster 
method, which is not rejection-free, depends on a tun- 
able parameter, n. In order to maximize the acceptance 
of the SW attempt (SWA), k should be chosen as close as 
possible to /3{e}(e). After every e change, we performed 
cycles consisting of 10^ Metropolis steps, n refreshing, 
then 10^ SWA, and a new k refreshing. The cycling was 
stopped, and k, fixed, when the SWA acceptance exceeded 
60%. Afterwards we performed 2 — 4 x 10^ SWA, taking 
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measurements every 2 SWA. In addition, we performed 
thermalization checks that included comparisons of hot 
and cold starts or even mixed configurations 

Our results for the latent- heat, Ae = — Cq, and the 
surface tension are in Fig. [31 The apparent location of the 
tricritical point (i.e. the p where both Ae and S vanish) 
shifts to upper p for growing L rather fast. For lattice 
sizes comparable with those of previous work, L — 16, we 
obtain Pt^^^ ~ 0.75, at a sizeable distance from p=l, 
but the estimate of pt increases very fast with L. 



all p < Pt, one expects that both ^{ed)/L and ^{eo)/L 
tend to non- vanishing and different limits for large Lj25j. 
On the other hand, for p > pt, 5(ed)/i is of order 
while ^(eo)/L ^ L^^^. For a fixed L, upon increasing 
p, the behavior goes from second-order like to first-order 
(see Fig[T]). Hence, a Finite-Size Scaling approach [l^ is 
needed. 
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FIG. 3: (color online). Top: Latent heat as obtained from 
/3(e) vs. spins concentration for several lattice sizes (lines are 
linear interpolations) . Data for p = 1 and L = 128 were taken 
from Ref. p^ . To illustrate the sample dispersion, we plot as 
well the scatter-plot of (N/L^ , A^^ye) for the 128 samples at 
L = 16 p = 0.85 and L = 64 p = 0.92. Bottom: as top part, for 
the surface tension. 

The PDFs for Ae and S, Fig. [21 display an interesting 
L evolution. When the /3(e) changes behavior from non- 
monotonic (L — 64, Fig. [Tl bottom-right) to monotonic 
(L = 128, Fig. [Tlbottom-left), the two PDFs becomes 
enormously wide|24|. see top panels in Fig[21 This arises 
because for many L = 128 samples, the curve /3{j}(e) 
is becoming fiat, or even monotonically decreasing (i.e. 
Ae = IJ = 0), while no such behavior was seen for L = 64. 
Only for p = 0.98, the width of the PDFs for Ae scales as 
i^^/^, as expected for a self-averaging quantity. Fig. [2]- 
bottom-left. The surface-tension is not self- averaging. 
Fig. [21-bottom-right. 

From Figs. [H [2] and [3] one cannot rule out that pt — I: 
a disordered first-order transition would not exist. For- 
tunately we can solve this dilemma by considering the 
correlation-length, obtained from the sample-averaged 
correlation function. 



FIG. 4: (color online) . Left: Correlation length in units of the 
lattice size, at phase-coexistence for the paramagnetic (top) 
and ordered (bottom) phases , as a function of spin concen- 
tration for several L (lines are cubic spline interpolations for 
data at fixed L). Right: Spin concentration where ^/L (data 
from left panel) coincide for lattices L and 2L versus 
see Eas. (|3l4p . Lines are a joint fit for x, pt, Ad and Ao. 

Consider the curves of ^(ed)/i versus p, for different L, 
Fig. m (left-top). There is a unique concentration, p^'^^ , 
where the correlation length in units of the lattice size 
coincides for lattices L and 2L. One has [261 



p^-^^^pt 



(3) 



An analogous result holds for ^{eo)/L. Since and 
are rather different, see Fig. [4l — right, a joint fit of all 
data yields an accurate estimate for the location of the 
tricritical point: 



Pt = 0.954(3), 



1.23(9), 
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dof 



4.23 



(4) 
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as ^2(g) ^ [-1 C(0,0,0)/C(27r/L,0,0)]/[2sin7r/L], 
where C is the Fourier transform of C(r) [rzllTst. 

We take the correlation-length in units of the lat- 
tice size at ed,eo as obtained from /3(e) (a jackknife 
method 18[ takes care of the statistical correlations) . For 



Of course, due to higher-order scaling corrections, Eq.Q 
should be used only for lattices larger than some 
^min -^^gj rpj^^ ^2 ^g^g acceptable taking Ljf'" = 16 
and L™™ — 12 (for the sake of clarity we do not dis- 
play data for L = 12 in the figures). We thus conclude 
that p = 0.98 is definitively in the first-order part of the 
critical line. 

We now look at ^/L at p^-^^, Fig[4l Consider C(ed)/i 
(f(eo)/i) as a function of {L,p), in the region p < pt- 
The salient features are: (i) for fixed L, f(ed)/-/j is a de- 
creasing function of p (^(eo)/i is increasing); (ii) for fixed 
p, ^(ed)/-L has a minimum {^{eo)/L has a maximum), at 
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a crossover length scale, Lco(p), that separates the first- 
order like behavior from the second order one; (iii) at 
the crossing point p^'^^ we have L < Lco{p^'^^) < 2L; 
(iv) at least within the range of our simulations, Lco(p) 
is a growing function of p. A standard scaling argument, 
combined with (i) — (iv), yields that ^{ed)/L at p^'"^^ is of 
order 1/ico {£,{eo)/L LcJ"^). If icob) diverges at pt, 
^{ed)/L at should tend to zero for large L, which 

is indeed consistent with our data. 

In this work, we have performed for the first time a de- 
tailed study of a disordered first-order transition in D = 3, 
by site-diluting the (5 = 4 Potts model, a system suffer- 
ing a prototypically strong first-order transition. A fairly 
small degree of dilution smooths the transition to the 
point of becoming second order, at a tricritical point, pt- 
A delicate Finite-Size Scaling analysis is needed to firmly 
conclude that pt < 1- We thus claim that (quenched) dis- 
ordered first-order transitions do exist in £) = 3, although 
quenched disorder is astonishing effective in smoothing 
the transition (we speculate that the percolative mecha- 
nism for colossal magnetoresistance proposed in [H could 
be fairly common in D — 3). We also observe that, for 
a given p < pt, a, crossover length scale Lco{p) exists 
such that for L < Lco(p) the behavior is first order like. 
The asymptotic second-order behavior appears only for 
L > Lcoip)- Our data are consistent with a divergence 
of Lco{p) at Pt. The successful location of the tricriti- 
cal point has been made possible by new definitions of 
the quenched average that avoids long-tailed PDF Q . It 
was crucial in this approach a recently introduced micro- 
canonical Monte Carlo method that features the entropy 
density rather than the free energy [l3| . 
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